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>v Abstract 



We define a new multispecies model of Calogero type in D dimensions with 
harmonic, two-body and three-body interactions. Using the underlying conformal 
SU(1, 1) algebra, we indicate how to find the complete set of the states in Bargmann- 
Fock space. There are towers of states, with equidistant energy spectra in each 
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tower. We explicitely construct all polynomial eigenstates, namely the center-of- 

O 

mass states and global dilatation modes, and find their corresponding eigenenergies. 
Qh! We also construct ladder operators for these global collective states. Analysing 



corresponding Fock space, we detect the universal critical point at which the model 



exhibits singular behavior. The above results are universal for all systems with 
Ctf ! underlying conformal SU(1, 1) symmetry. 
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1 Introduction 



The (rational) Calogero model describes N identical particles on the line which inter- 
act through an inverse-square two-body interaction and are subjected to a common 
confining harmonic force. Starting from the inception [1], the model and its vari- 
ous descendants (also known as Calogero-Sutherland-Moser systems [2]) continue to 
be of interest for both physics and mathematics community, primarly because they 
are connected with a number of mathematical and physical problems, ranging from 
random matrices and symmetric polynomials [3] to condensed matter systems and 
black hole physics [4]. The model is also connected to Haldane's exclusion statistics 
[5]. The role of Haldane statistical parameter is played by (universal) coupling con- 
stant in the two-body interaction. In Haldane's formulation there is the possibility 
of having particles of different species with a mutual statistical coupling parame- 
ter depending on the species coupled. This suggest the possible generalization of 
the ordinary one-dimensional Calogero model with identical particles to the one- 
dimensional Calogero model with non-identical particles. Distinguishabillity can be 
introduced by allowing particles to have different masses and different couplings to 
each other. In this way a one-dimensional multispecies Calogero model is obtained 
[6,7]. 

Further generalization can be achieved by formulating the model in dimensions 
higher than one. In the case of single-species model(s), it was shown that some 
exact eigenstates (including the ground state) can be obtained for a D dimensions 
provided that a long-range three-body interaction is added [8]. The inevitable ap- 
pearance of three-body interaction in D > 1 makes any analysis of such a model (s) 
highly nontrivial and very little is known about their exact solvability. Some progress 
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has been achieved only recently for a class of two-dimensional models with identical 
particles [9]. 

In a present Letter we propose a new type of partially solvable multispecies model 
of Calogero type in D dimensions. In addition to the harmonic potential, it con- 
tains two-body and three-body interactions with coupling constants depending on 
the particle's species. We also allow particles to have different masses. In this 
way we incorporate both generalizations mentioned above into a single model. We 
indicate how to obtain (in principle) all eigenstates of the model Hamiltonian in 
Bargmann representation. The spectrum of states shows a remarkable simplicity. 
There are towers of states with equidistant energies. We are able to find all poly- 
nomial eigenstates and corresponding eigenenergies of the Hamiltonian, describing 
global collective states. Closer inspection of the Fock space, corresponding to the 
relative motion of particles, reveals the existence of the universal critical point at 
which system exhibits singular behaviour. This result generalizes that mentioned 
in [7,10]. Our results are universal and applicable to all systems with underlying 
£77(1, 1) algebra. 

2 A model Hamiltonian 

We start the analysis with observation that the exact wave functions of the Calogero 
model are highly correlated. These correlations are encoded in the wave functions in 
the form of a Jastrow prefactor (x« — Xj) v for any pair of particles The exponent 
of the correlator is related to the strength of the two-body interaction. It is then 
plausible to make an ansatz for the most general ground state wave function for the 
N distinguishable Calogero-like particles in D dimensions in the form (h = 1) 
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M/ (ri,...,rV) =Ae-^L^ (1) 
where the Jastrow prefactor is generalized to 



A = n|ri-^p, ^ = ^, i,j = l,2,...iV. (2) 

Here, mj are masses of the particles, is the frequency of the harmonic potential and 
Uij are the statistical parameters between particles i and j. In principle, one could 
start with any wave function with no nodes, except at the coincidence points, and 
which is continuosly connected with Gausss function when parameters — > 0. Note 
that for — u, rrii — m and D ^ 1, Eq.(l) smoothly goes to exact ground state 
of the Calogero-Marchioro model [11], so the wave function (1) is a natural choice. 
Adopting the reasoning from Ref.[12], we can ask for what kind of Hamiltonian is 
the wave function (1) the exact ground state. It turns out that ^o(rt, ■••,r/v) will 
be, for sufficiently small deformations u^, the exact ground state of the Hamiltonian 

1 _ v ij (v ij + D-2) 1 J_ 
1 VijVik(ri - rj)(r; - f k ) 

2 i#i,i#fc \n-r k \ 

such that 

tf^o = Wo, (4) 

E o = + ^ ftf) = ^e - (5) 

The ground state (1) and the Hamiltonian (3) are invariant under the group of per- 
mutation of N elements, Sat, generated by exchange operators Kij [13]. Operators 



Kij interchange indices % <-> j in all quantities, i.e. rrii <-> m^, ^ <-> ^jfc, fj <-> 
Pi <-> Pi- 

For £) = 1 the three-body term in (3) identically vanish if z/y = i/, rrii = m or if 
v,ij = cmiirrij, a beeing some universal constant [7]. Unlike in one dimension, how- 
ever, it does not vanish in higher dimensions and plays a crucial role in the analysis 
that is to follow. 

Let us perform the non-unitary transformation on ty , namely — A _1 ^ r - It 
generates a similarity transformation which leads to another Sn invariant Hamilto- 
nian H = A' 1 HA. We find H as 



= u 2 T + -T_. (6) 

where we have introduced 

i=l 1 i<j \ r i~ r j\ " tl "h 

Y n y N 

T+ = t £ wf, T o = -(E^V, + e ). 

z i=i z i=i 

The operators T±, T satisfy the ST/^l, 1) algebra 

[T_,T + ] = 2T , [T ,T ± ] = ±T ± . (8) 
The following identity (i.e. similarity transformation) holds for ui ^ 0: 

# = uj 2 T + - T_ = 2u;,S T S' 1 , (9) 
5 = e- T +e-^ T -. 
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Owing to this identity, we can employ Bargmann representation and construct iter- 
atively Bargmann- Fock space of eigenstates. We begin with state $ , which is the 
lowest weight vector of the operator T_ and also an eigenstate of T : 

T_$ = 0, T $ = |$ . (10) 

In our case, $o = 1 and eo is given in Eq.(5). 

The tower of excited states ( level 0-tower ) is obtained by succesive application of 
T_ operator: 

T_$ 2p = $ 2p _ 2 , 2uT^ 2p = u{2p+e )<5> 2p , p = 0, 1, 2, 3.... (11) 

The states $ 2p are either polynomials or irrational functions of homogenity 2p, and 
are eigenstates of T . Two succesive states differ in energy by an amount 2uj. 
Similarly, one can construct towers of states at level 1, $ 2 p +i , p — 0, 1, 2..., 
h = (h = 1, 2, ...N; ai = 1, 2, ...D), using 

T_<^ = 0, T_^ p+1 = 2uT ^ p+1 = u(2p + e[^i p+1 . (12) 

Here, ef 1 is energy of the first excited state which tends to (1 + 4^) in the limit 

Uij — > 0. Two succesive states also differ in energy by an amount 2u. 

Following the procedure, one gets the towers of states at level k, < k < ND, using 

T-*t -' h = 0, T. = $W- 2) 2uT ^ 2p ^ = o;(2p+^'- A )$fe A . 

Here, the energies e I k 1, "' ,Ik tends to (fc + 4p) in the limit ^ — > 0. The states 
< ^ ) 2p+fc 7fe are eigenstates of the Hamiltonian 2uT = S^ 1 H S, Eq.(9). Particularly, 
the state $ = 1 is a ground state (i.e. the lowest energy eigenstate) for all towers 
if e < e I k 1, "' Ik , VJi, ...,/fc and for all indices fc. 

Notice that the operator T + of Eq.(7), acting on the particular state in the given 



tower, gives an another state in the same tower with energy greater by an amount 
2uj (see also Sec. 3, Eq.(17)). 

The procedure outlined above in Eqs. (10-12) is exaustive, i.e. it gives all eigenstates 
of the S'-transformed Hamiltonian S^ 1 H S (cf. Eq.(9)) , provided one is able to 
solve differential equations (11) and (12). This is a non-trivial task, even in D = 1. 
However, one can readily show that this procedure, when applied to the system of 
N D-dimensional free harmonic oscillators (i/^ = in Eqs. (6) and (3)) yields the 
following set of eigenstates for H = H: 

N D 

S'dilKr), 7^ = 0,1,2...,, (is) 

i=l a=l 

where (cf. Eqs. (7) and (9)) 

1 x ^JV ,2 1 1 y^JV l 

S = e + e"^ = e _tJ 2 z^=i m » r ; e s ^ i=1 ~ 1 . 

The corresponding eigenenergies are 

ND N D 

"(-5- + EE **.«)■ ( 14 ) 

z i=l a =l 

For the convenience of the reader, we describe in the next Table the structure of the 
few lowest towers of states at level k (Eqs.(10)-(12)) in this simple case. 



k 


level- A; tower 


indices 





1 




1 


n 


/ = (i, a) 


2 


rhri 2 




3 


r h r l2 r h 











For the general case, Eqs.(10)-(12), the towers of states at level k, S^ 1 ''"' 7 *, need not 
have such a simple monomial structure since they can be , in principle, homoge- 
nious irrational functions. From Eq.(13) one can count the number of states at each 
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level of given homogenity. For example, there are ND states of homogenity one, 
ND + — ^- p ~ 1 ) states of homogenity two, etc. There are 2 nd towers in total. 
Now, one can put an interesting question, namely is there, for sufficiently small 
deformation parameters Vij, one-to-one correspondence between our multispecies 
model H(i/ij), Eq.(6), and iV D-dimensional free harmonic oscillators 
H(uij = 0) = H{vij = 0). According to our analysis, there is no unique similarity 
transformation between these two systems. However, there is similarity transforma- 
tion between given tower in the interacting system [y^ ^ 0) and analogous tower 
in the free system {y^ = 0), up to the constant e I k 1 '"' Ik . Particularly, this was shown 
for D — 1 and identical particles (z^ = v) in Ref. [9]. In that case, the eigenstates 
are restricted to Sjv-symmetric representations. 

We are unable to find towers of states by solving differential equations (11-12) 
in general. However, as we will show in the next section, we are able to construct 
global collective states for the Hamiltonian (6). These states represent all states of 
the polynomial type in Bargmann representation in generic case. Moreover, these 
states are universal for all systems with underlying conformal SU(1, 1) symmetry. 

3 Ladder operators and Fock space representa- 
tion for global collective states 

It is convenient to introduce the center-of-mass coordinate R and the relative coor- 
dinates pi [14]: 

Y n N 

m i=i i=i 

— * -* — * Tfl ' — * 

ik = ?i-R, V ft = V i - T pV Jl (15) 
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They satisfy identities X^Li m>iPi = Z)£Li V Pi = 0. In terms of the variables just 
introduced, the Hamiltonian H and wave function separate into parts which 
describe center-of-mass motion (CM) and relative motion (R), namely H = Hcm + 
H R and *o(n,...,nv) = *oCR)*o(Pi,...,Pat). 

Using Eqs.(7) and (15) we define creation (+) and annihilation (-) operators 

A, ± = h— +u,T+)tT , (16) 

Z Ul 

which satisfy the following commutation relations (a, (3 = 1, 2, • • ■ D) : 

[AC, AC] = AC, [A 2 ,A 1 + ] = AC, 

[A 2 ,A 2 + ] = -, [H,A 1 ± ] = ±uA 1 ± , 

[H,A 2 ± ] = ±2luA 2 ± . (17) 

Notice that A 2 ± = ST±S~\ AC = SRS' 1 and AC = SVrS' 1 , with S defind 
in Eq.(9). They act on the Fock vacuum |0) oc S&o(ri, ...,rW) as 

AC\0)=A 2 \0) = 0, (0|0) = 1. (18) 

The excited states in the Fock space, corresponding to global collective states, are 
of the form 

Atc i ' i ---A+ D ni > D Ar\o) = n w ni ^rio), m 

a=l 

where n^ a — 0, 1, 2...(Va) and n 2 — 0, 1, 2... 

The repeated action of the operators on the vacuum |0) reproduces, in the coor- 
dinate representation, Hermite polynomials H ni a (R a \J Muj). Similarly, the repeated 



action of the operator A\ on the vacuum |0) reproduces hypergeometric function, 
which reduces to associated Laguerre polynomials L^°^ _ 1 (2o;T + ) for certain values 
of parameters. The states (19) are eigenstates of the H with the energy eigenvalues 
(cf. last two equations in Eqs.(17)) 

D 

E nia . n2 =uj{^n^ a + 2n2 + SQ). (20) 

a=l 

This is the part of the complete spectrum which corresponds to center-of-mass states 
and global dilatation states, respectively. 

Now we show that the states (19) are perfectly normalizable,i.e. quadratically 
integrable and physically acceptable for both Hamiltonians H and H, provided that 
e > y. First, we completely decouple CM- and R- motion by introducing another 
set of the creation and annihilation operators {B^B^}: 

Bf^A^- 1 -^, (21) 

such that 

[Ax^fla*] =0. (22) 

Hence, we get 

H R = uj[B 2 , B+], [H R , B±] = ±2cuB±, 

HcM = lu;J2{A^,A ha + } + . (23) 

Z Q=l 

The Fock space now splits into the CM-Fock space, spanned by Y\a=i{ J ^t,a) ni ' a |0)cm 
and the R-Fock space, spanned by £>^" n2 |0)R, where \0)cm oc e~2 MR2 and \0)r oc 
e ^' ra,,, «. We point out that the R-modes are universal for all systems with 
underlying conformal SU(1, 1) symmetry, i.e. for the Hamiltonians of the form 
H = -T_ + u 2 T + + 7 T , where T±,T satisfy 5*7(1, 1) algebra (8). 
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Closer inspection of the R-Fock space of the Hamiltonian Hr, Eq.(23), reveals the 
existence of the universal critical point defined by the zero-energy condition 

^=(^ + ^ = 0. (24) 

At the critical point the system described by Hr collapses completely. This means 
that the relative coordinates, the relative momenta and the relative energy are all 
zero at this critical point. There survives only one oscillator, describing the motion 
of the centre-of-mass. Such behaviour resembles some features of the Bose-Einstein 
condensate. It was first noticed in Ref.[9] for the case D — 1, v i3 ■, — v and = m. 
In that case the critical point (24) is simply at v — (Notice that there is also 
critical point at v = 1 + for this case). Of course, for the initial Hamiltonian H, 
Eq.(3), which is not unitary (i.e. physically) equivalent to H, this corresponds to 
some Vij < 0, satsfying Eq.(24), and the norm of the wave function (1) blows up at 
the critical point. For z^- negative but greater than the critical values (24), the wave 
function is singular at coincidence points but still quadratically integrable. Out 
of the critical point we have one-to-one correspondence between our multispecies 
system (6) and the system of N D-dimensional free oscillators, at least for the 
dilatation states B^ U2 \0) R . 

4 Conclusion 

In summary, we have defined a nontrivial many-body Hamiltonian H (Eq.(3)) of 
Calogero type in D dimensions with two- and three-body interactions among non- 
identical particles. Strength of the interactions, i/^, depends on the particle's species 
and this feature makes any analysis of such a model nontrivial, even in D = 1. Us- 
ing underlying SU(1, 1) structure of the transformed Hamiltonian H (Eq.(6)) and 
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Bargmann representation we outlined a procedure which gave in principle all eigen- 
states of the Hamiltonian. While we were unable to solve corresponding differential 
equations (11,12), we were able to find some general features of the solutions. There 
are towers of states with equidistant energy spectra. In each tower two neighbour- 
ing states differ in energy by 2oj. Moreover, we managed to solve H partially, i.e. 
we explicitely found its global collective states, corresponding to the center-of-mass 
motion and the relative motion of particles. Those are all polynomial solution in the 
Bargmann representation in generic case. We also found their eigenenergies. The 
spectrum of collective modes, Eq.(20), is linear, equidistant and degenerate. It is 
also found that, for u ij — — (N — 1)D, the Fock space, corresponding to the 
relative motion of particles, contained states of zero norm and the whole system 
exhibited singular behaviour. At this critical point the ground state wave function 
of the Hamiltonian H, Eq.(3), posseses infinite norm. 

If we consider identical Bose (Fermi) particles, with = m, and = u, the 
eigenstates are restricted to SAT-symmetric (antisymmetric) functions and the criti- 
cal point is at v = — -j^. Our analysis of multispecies Calogero model gives deeper 
insight on the single-species Calogero models in higher dimensions. 
All results presented here are common and universal for all systems with underlying 
conformal SU(1, 1) symmetry. The potentially most interesting applications of our 
results might be in two dimensions and quantum Hall efect. 
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